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ABSTRACT 

Studies of some bright, super-Eddington transient pulsars show a negative correlation between the energy of the cyclotron resonance 
scattering feature (CRSF) and the bolometric luminosity. For Her X-1, using repeated RXTE observations during 1996-2005, the 
inverse dependence was found: the energy of the cyclotron line increases as the luminosity increases. The X-ray flux measured by 
the RXTE/ ASM (2-10 keV) has been assumed to represent the luminosity - more precisely: the maximum X-ray flux reached during 
the respective 35 d Main-On. Here, we question whether the ASM flux is really an accurate measure of the bolometric luminosity 
of the source. We redetermined the energy of the cyclotron line and performed spectral fits using the combined data from the PCA 
(3.5-60 keV) and HEXTE (20-75 keV) instruments on RXTE of the same 35 d cycles as used in the original work to determine the 
bolometric flux from those spectra. We confirm the result of the original analysis that the cyclotron line energy changes by ~7% for a 
change in flux by a factor of two. 
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1. Introduction 

Her X-1 is one of the most observed and well-studied accreting 
binary X-ray pulsars. Since the first observation made by Uhuru, 
the 35 d periodicity of the source is well known. The X-ray light 
curve shows two o«-states (high X-ray flux) and two Oj^-states 
(low X-ray flux), with a Main-On {~ 7 orbital periods) and a 
Short-On (~ 5 orbital cycles), separated by two Oj^-states (~ 4-h5 
orbital cycles each). The maximum X-ray flux of the Main-on is 
higher than the maximum X-ray flux of the Short-On by a factor 
of three to four. 

The modulation of the X-ray flux is believed to be caused 
by the periodic obscuration of the X-ray source by the pre- 
cessing disk, which is b elieved to be inclined and warped 
dGerend & Bovntonll976h . The onset of the flux, the turn-on (of- 
ten identified with 35 d phase 0.0), is believed to occur when the 
outer rim of the disk opens up the view to the X-ray emitting re- 
gions near the polar caps on the surface of the neutron star, while 
the flux decrease towards the end of the o«-states is associated 
with the inner parts of the disk covering these regions from the 
observer. 

Her X- 1 is also the first X-ray pulsar for which a cyclotron 
line was discovered dTriimper et al.ll 19781) . The discovery of this 
line-like feature played a key-role in the measurement of the 
magnetic fields, providing the first direct measurement of the B- 
field of a neutron star. The energy of the cyclotron line is related 
to the magnetic field by the formula B12 - (1+z) Ecyc/11.6keV, 
where B12 is the magnetic field strength in units of 10'^ gauss, z 
is the gravitational redshift and Ecyc is the energy of the cyclotron 
line. This feature is now referred to as a cyclotron resonant scat- 
tering feature (CR SF) and seems to be quite common in accret- 
ing X-ray pulsars dCoburn et al.ll2002h . The cyclotron line is an 
absorption feature produced by the resonant scattering of pho- 
tons on electrons. In the ~ 10'~ gauss magnetic field, the elec- 
trons are in quantized energy states (with respect to their move- 
ment perpendicular to the magnetic field), the so-called Landau 
levels. The energy levels are nearly equidistantly spaced and 



photons with energies equal to n times the fundamental Landau 
energy could take part in this scattering. 

In a few transient pulsars such as V0332+53 and 4U 
01 15+63, a negative correlation between the CRSF and the bolo- 
metric luminosity of the source has been observed: the energy of 
the cyclotron Une decreases as the X-ray luminosity increases 
(Miharaetal. 1998:M owlavi et alJl2006tlNakaiima et al.ll2006l: 
Tsvgankov et al. 2006) . Her X-1, however, shows the opposite 
behavior. Repeated measurements of the cyclotron line energy 
with different instruments such as RXTE and INTEGRAL in 
the X-ray spectrum of Her X-1 have revealed a positive correla- 
tion between the (pulse phase-averaged) cycl otron line energy 
Ecyc and the X-ray flux dStaubert et al.l uOOTl) . more precisely 
with the maximum X-ray Eux observed during the Main-On of 
the respective 35 d cycle, as measured by RXTE/ ASM in the 2- 
lOkeV range, which is assumed to be characteristic of the cur- 
rent accretion state and luminosity of the source. It has, however, 
been questioned whether the 2-lOkeV flux can really be taken 
as representative of the bolometric luminosity. 

In this paper, we attempt to answer this question by re- 
analysing the RXTE observa tions of the same 35 d cycles as 
used in the original analysis of lStaubert et alJ d2007h . Using data 
from both instruments (PCA and HEXTE), we re-determine the 
cyclotron line energy and measure the 3.5-60 keV bolometric 
X-ray flux by performing spectral analyses of eleven Main-Ons 
between 1996 and 2005. We show that the 2-lOkeV flux is a 
good measure of the bolometric flux and that the positive cor- 
relation between the phase-averaged cyclotron line energy Ecyc 
and the maximum 35d Hux, or the luminosity, is confirmed. 

2. Observations 

Her X-1 has been repeatedly observed by RXTE since 1996. We 
analyzed observations of eleven 35 d Main-Ons for which there 
was photon statistics of sufficient high quality to allow a spectral 
analysis. Those cycles were those of numbers (nos.) 257, 269, 
303, 304, 307, 308, 320, 323, 324, 343 and 351. For the defini- 
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Table 1. Details of RXTE observations of Her X-1 used for the 
spectral analysis. 



Observation 
month/year 



35 d Main-On 
cycle number ' 



Center 
MJD 



July 96 


257 


50029.75 


September 97 


269 


50707.06 


December 00 


303 


51897.69 


January 01 


304 


51933.67 


May 01 


307 


52035.48 


June 01 


308 


52071.16 


August 02 


320 


52492.96 


November 02 


323 


52599.32 


December 02 


324 


52634.01 


October 04 


343 


53300.95 


July 05 


351 


53577.35 








' Cycle numbering 


accordins tolStaubert et alj (120091) 



tion of cycle counting we refer to IStaubert et al.l (12009*) (cycle 
no. 313 was not included, despite the high quality statistics, be- 
cause no ASM observations were available). We cover a period 
often years of observation from July 1996 (no. 257) to July 2005 
(no. 351), see Table [1] Here we used data from both RXTE in- 
struments: PC A in the energy range 3.5-60keV, and HEXTE in 
the energy range 20-75 ke V. 



3. Spectral analysis 

For the spectral analysis, we used XSPECQ (12.6.0) and gener- 
ated spectra for all observed Main-Ons, summing all available 
data. For each spectrum (one for each Main-On), two quantities 
were determined: the central energy Ecyc of the cyclotron ab- 
sorption feature and the integrated flux in the range 3.5-60keV 
(in units of erg cm"^ s"'). As the spectral function, we used the 
highecut' model (based on a power-law continuum with an ex- 
ponential cut-off) and a multiplicative Gaussian absorption line 
for the cycl otron resonan t scatte ring feature (CRSF), in the same 
way as e.g. lCoburn et al.l(l2002h and lStaubert et al.l (l2007h . Data 
from both RXTE instruments were used: PCA (PCU2 only) in 
the energy range 3.5-60keV, and HEXTE in the energy range 
20-75 ke V. Cold material absorption was taken into account in 
each fit. All values were consistent to within two standard de- 
viations with the mean value of Nh - 1.08 x 10^^ cm"-. We 
note that in the original analysis PCA data were used only up 
to 25 keV to define the continuum at the lower energies. New 
response matrices are now availablq^ that allow us to use the 
PCA up to ~50 keV, but we found that the PCA can indeed be 
used up to 60 keV, thereby contributing information about the 
cyclotron line around 40keV. This choice is confirmed by the 
agreement between the PCA and HEXTE spectra in the over- 
lapping region of the two ins truments (40-60 keV). T his is also 
confirmed in the analysis by iRothschild et al.l (1201 ll) of RXTE 
observations of Cen A in which PCA data were successfully 
used up to 60 keV. In the new spectral fits, we also added sys- 
tematic uncertainties of 0.5% (which are recommended for use 
wi th the new r esponse matrices^), while in the original analysis 
bv lStaubert et al. (2007) 1% uncertainty was added. 




'++% M'''^^'''tT#|||r|i|uk 



' http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec 

^ http://www.universe.nasa.gov/xrays/programs/rxte/pca/doc/rmf/ 



pcarmf-11.7 



Energy (keV) 

Fig. 1. Example of a spectral fit to an RXTE observation of 
Her X-1 in July 1996 (cycle no. 257). The observation is cen- 
tered at MJD 52599.36 and the integration time is 31 ks for PCA 
and lOks for each of the two HEXTE clusters. Black: PCA, red: 
HEXTE-A, and green: HEXTE-B, respectively; top: count rate 
spectra; middle: residuals with respect to a fit of a continuum 
model; bottom: residuals with respect to a fit which includes a 
cyclotron line. 



For this spectral re-analysis, only data of RXTE PCU2 were 
used to ensure maximum uniformity in the data set, analyzed 
in a uniform way. Fig. 1 shows an example of a spectral fit 
to data of July 1996 (cycle no. 257). The observation is cen- 
tered at MJD 52599.36 and the total integration time is 31 ks for 
PCA/PCU2 and lOks for each of the two HEXTE clusters. 

The bolometric flux for the individual spectra was found by 
integrating the fit function (over the 3.5-60 keV range), using 
the XSPEC routine flux'. The X-ray flux of Her X-1 varies as 
a function of phase of the 35 d modulation because of variable 
absorption by the accretion disk. 

The sampling of the various Main-Ons in groups of pointed 
observations was rather different, and the fluxes determined from 
the spectra are mean values of the observations, which happened 
to be made at different phases of the 35 d modulation. These 
fluxes cannot directly be compared to the maximum X-ray flux 
of the respective Main-On. To find the comparable bolometric 
fluxes in the 3.5-60keV range that represent the maximum flux 
for the particular cycle, it was necessary to scale the fluxes found 
in the individual spectra to the maximum flux. The scaling factor 
is the ratio of the maximum flux to the mean flux of the complete 
light curve, that is of the same data used to generate the spectra. 
To find the maximum flux, all observed RXTE/PCA light curves 
were fitted by a function representing a template of the mean 
35 d Main-On modulation, the shape of which was taken from 
an overl ay of many 35 d light curves observed by RXTE/ ASM 
(see e.g. lKlochkov et al1l2006l) . For this procedure, we are only 
interested in the overall 35 d modulation, that is without eclipses 
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37 38 39 40 41 42 

Cyclotron line energy from Stoubert et al. 2007 (keV) 



4 5 6 7 8 

ASM max. flux from Stoubert et ol. 2007 (cts/s) 



Fig. 2. hefv. Ecyc from the spectral re-analysis versus Ecyc from the original analysis (IStaubert et al.ll2007h . Right: Scaled maximum 
3.5-60keV flux (in units of 10"^ erg cm"^ s ') from the spectral re-analysis versus the maximum 2-lOkeV ASM flux (in units of 
cts s ') from the original analysis. The Pearson correlation coefficients are 0.95 (left) and 0.82 (right). 

Table 2. Phase-averaged cyclotron line energies Ecyc and maximum fluxes reproduced from lStaubert et all (l2007h (1) and from the 
spectral re-analysis of this work (2). 



35 d cycle 
number' 



257 
269 
303 
304 
307 
308 
320 
323 
324 
343 
351 



E™(1) [keV] 



fStaubert et al. (2007^) 



max. flux (2-10 keV) 

ASM [cts s-'] 
fStaubert et al. f2007')') 



E,y,(2) [keV] 
(this work) 



max. flux (3.5-60 keV)- 

[10"' erg cm"- s"'] 

(this work) 



41.12: 

40.62 : 
40.07 : 
39.05 : 
39.93 : 
39.73 : 
40.01 : 

40.51: 
40.60 : 
38.51 : 

38.95 : 



0.55 
0.49 
0.31 
0.55 
0.63 
0.52 
0.29 
0.13 
0.41 
0.51 
0.52 



7.37 ± 0.34 

7.49 ± 0.73 
6.04 ± 0.47 
5.72 ± 0.34 
7.15 ±0.50 
6.93 ± 0.20 
7.19 ±0.26 
7.64 ± 0.30 
7.55 ± 0.34 

4.50 ± 0.24 
5.12 ±0.37 



40.61 : 
40.45 : 
39.81 : 
38.72: 
39.78 : 
40.04 : 
39.94 : 

40.16: 

40.37 : 
38.26: 
39.28 : 



0.34 
0.41 
0.24 
0.31 
0.49 
0.33 
0.29 
0.11 
0.35 
0.43 
0.60 



7.99: 
8.02: 
7.27: 
7.14: 
8.26: 
8.18: 
9.55: 
9.37: 
9.16: 
7.05: 
6.61: 



0.45 
0.46 
0.41 
0.40 
0.47 
0.46 
0.36 
0.44 
0.52 
0.40 
0.38 



' Note that for 35 d cycle numbers 303 and larger the corresponding numbers in lStaubert et alj ( 120071) are larger by 1 . This reflects the observation 
that there must have been an extra cycle during the long anomalous low before cycle 303. However, using the numbers as given here allows to use 
them in an ephe meris f or a rough prediction of the 35 d turn-ons using a mean period of 35.88 d. We do not doubt the physical reality of the extra 
cycle found by _Staubert_et al. (2007). 
^ The flux in Col. 5 represents the 35 d maximum flux in the 3.5-60 keV range (see text). 



and dips, which was described by the analytical function: 

1 



/(0=Aol- 



1 +exp(Bi(/-Ci))/ 1 +exp(B2(f-C2))' 



under the condition that negative functional values are set to 
zero. Here /(f) is the flux as a function of time f in MJD, Aq is the 
ampJifude in units of (3-20keV) PCA cts s"' per PCU,Ai = 1.25, 
Bi = 1.87, and ^2=0. 80 are fixed functional parameters, Ci is the 
time reference in MJD, and C2=Ci -1-8.4 d (this offset represents 
the length of the standard Main-On). The time of reference Ci 
was fixed to the MJD of the turn-on that had been determined 
from the corresponding ASM light curve. The amplitude Aq is 
the only free-fitting parameter. The maximum of the fitting func- 
tion (on average found to be equal to Ao-5.4 cts s"') is then taken 
as the maximum Hux (in units of PCA cts s"') for this cycle. 
Dividing this maximum flux by the corresponding mean flux of 
the complete light curve (containing all the photons that are also 
used in generating the corresponding spectrum) leads to a scal- 
ing factor. Multiplying the 3.5-60 keV flux found from the spec- 



tra by this scaling factor then leads to the maximum bolometric 
flux (in units of erg cm~^ s"'). These bolometric flux values can 
now be compar ed to the maximum A SM flux values used in the 
earlier work bv lStaubert et al.1 (l2007h . 

To test whether the above procedure in determining the 
maximum bolometric fluxes could be influenced by systematic 
changes of the spectral shape over the 35 d phase, we have inves- 
tigated spectra of cycle no. 323 which provides a good coverage 
of a complete Main-On. Systematic changes were only found in 
one spectral parameter, the power-law index, with a magnitude 
of +3%, which corresponds to a negligible relative change of the 
integrated flux of less than 5 x 10"^%. 

The uncertainties in the bolometric flux values are dominated 
by the uncertainties in the scaling factors, which in turn are dom- 
inated by systematic uncertainties in determining the maximum 
Main-On flux in fitting the observed light curve by the common 
analytical function. These systematic uncertainties are estimated 
to be 3%, the main contributions being the uncertainty in the 
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Fig. 3. hefv. Ecyc from the spectral re-analysis versus the scaled maximum 3.5-60keV flux (in units of 10 erg cm s ). Right: 

^cyc 



Ecvc from the spectral re-analysis versus the maximum 2-lOkeV ASM flux (in units of cts s ') from the original analysis. 



turn-on time and the fluctuations in the shape of the modulation 
from cycle to cycle. 

Table 2 compares t he results of this re-a nalysis to those of 
the original analysis bv lStaubert et alJ (l2007h . 



4. Results and discussion 

In Fig.|2](left), we plot the newly determined cyclotron line en- 
ergies against those from the original analysis. Overall, we find a 
very good agreement between the new values for Ecyc and those 
from the original analysis; considering the respective uncertain- 
ties, the differences are smaller than 0.79 standard deviations for 
ten of the eleven values and 2. 1 standard deviations for one value 
(cycle no. 323). 

In Fig.|2](right), the 3.5-60keV maximum bolometric fluxes 
are plotted against the corresponding ASM fluxes of the origi- 
nal analysis. There is a good linear relationship between the two 
fluxes: the slope (taking the uncertainties in both variables into 
account) is 0.88+0.15 (10"'' erg cm^^ s"')/(ASM cts s"'), and 
the Pearson correlation coefficient is 0.82, corresponding to a 
chance probability of no correlation of P< 10"^^. This demon- 
strates that the maximum ASM flux can be taken as a good mea- 
sure of the maximum bolometric X-ray flux (and luminosity) of 
Her X-1 during the respective 35 d cycle. The variation in flux 
from one 35 d cycle to the next is small, such that the maximum 
observed flux of a particular 35 d Main-On can be considered a 
good measure of the luminosity of the source during this partic- 
ular cycle. This is why the maxim um ASM flux was use d as a 
reference in the original analysis bv lStaubert et al.l (l2007h . 

The final correlation between the cyclotron line energy and 
the X-ray flux is given in Fig. [3] in two ways: we correlate 



Ecyc from the re-analysis with the scaled 3.5-60keV maximum 
Main-On flux (Fig. [3] left), as well as the maximum ASM flux 
(Fig. |3] right). The corresponding slopes of the linear fits to 
these data (taking the uncertainties of both variables into ac- 
count) and the corresponding Pearson correlation coefficients r 
are: (i) for the Main-On flux presented in Fig. [3] (left) slope = 
0.62+0.19 (keV)/(10-^ erg cm^^ s"') and r = 0.62 (P=2 x lO'^); 
(ii) for the maximum ASM flux presented in Fig. [3] (right) slope 
= 0.67+0.14 (keV)/(ASM cts s"') and r = 0.90 (P^IO^*). 

The correlation seen in Fig. |3] (left) is somewhat less con- 
vincing than that of Fig. [3] (right) (and that in the original anal- 



ysis of IStaubert et al.l r2007h . We attribute this to the unavoid- 
ably larger uncertainties associated mainly with the scaling of 
the bolometric flux measured for the individual spectrum to the 
flux that does describe the maximum flux of the particular 35 d 
cycle. The originally used maximum ASM fluxes are, in con- 
trast, simple direct measurements. 



5. Summary 

We have re-analyzed observations of Her X-1 in its Main-On 
state by RXTE between 1996 and 2005 with respect to its X-ray 
spectrum. Using data from both instruments (PCA and HEXTE) 
we performed a spectral analysis of observations of eleven Main- 
Ons and determined the cyclotron line energy Ecyc and the 3.5- 
60 keV flux for each of those Main-Ons. This observed flux was 
then scaled to a flux representing the maximum flux of the par- 
ticular Main-On. We co nclude that the maxim um ASM flux used 
in the original analysis (IStaubert et al.ll2007n can really be taken 
as a measure of the luminosity of the source because it scales 
well with the maximum 5-60keV flux estimated through the 
spectral analysis. This is evident from Fig. |2] (right) that shows 
there is a good linear relationship between the bolometric 3.5- 
60keV flux with the 2-10 keV flux measured by RXTE/ ASM 
(both fluxes refer to the maximum 35 d Main-On flux). 

The information contained in Figs. 2 and 3 (left and right) 
provides an internally consistent picture. We consider the com- 
bined evidence of all the correlations shown as proof of the 
correctness of the positive correlation between the cyclotron 
line en ergy and source luminosity as suggested in lStaubert et al.l 
(l2007l) . A direct comparison to the original analysis can be made 
by considering Fig. |3](right): the slope of the linear best fit is de- 
termined to (0.67+0.14) keV/(cts s"'), which is in good agree- 
ment with the (0.66+0.10) keV/(cts s"') found in the original 
analysis. The final and confirmed result with respect to this cor- 
relation can then be stated as follows: the value of the cyclotron 
line energy Ecyc increases by ~7% for a change in flux of a factor 
of two. 

In conclusion, we briefly describe the physical significance 
of the observed dependence - of either sign - of the cyclotron 
line energy with source luminosity. The negative correlation 
(decrease of Ecyc with luminosity), which has been repeatedly 
observed for high luminosity transients such as V 0332+53 
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and 4U 0115+63 (iMihara et alJ [19981 iMowlavi et alJ l2006t 



iNakaiima et al.l l2006t iTsvgankov et al.l l2006h . has been inter- 



preted in the following way: when the mass accretion rate (and 
hence the luminosity) increases, the height of the radiative shock 
above the neutron star surface increases, leading to a decrease in 
the effective magnetic field strength in the scattering region and 
therefore to a decrease in Ecyc- This is in line with theoretical 
conside rations about accret ion in the (locally) super-Eddington 
regime ( iBurnard et al.ll99ll) . In the sub-Eddington regime, how- 
ever, believed to be realized in Her X- 1 , the deceleration of the 
accreted material is tho ught to be due to Co ulomb drag and col- 
lective plasma effects dNelson et al.lll993l) . In this case, under 
an increased accretion rate the atmosphere is compressed by the 
ram pressure of the infalling material and the scattering region 
moves closer to the neutron star surface. This is equivalent to an 
increase in effective field strength and an increase in Ecvk, hen ce 
to a positive correlation with luminosity dStaubert et al.ll20()7b . 

We would also like to add that the recent analysis of 
iKlochkov et al.l (l201lh supports the physical correlation between 
the cyclotron line energy and the luminosity of three sources: 
4U 0115+63, V 0332+53, and Her X-1. In this analysis data 
from a short time interval were used, and spectra were gener- 
ated by summing up photons belonging to individual pulses in 
selected ranges of pulse amplitude. In this pulse-to-pulse vari- 
ability study, the variations in X-ray flux (source luminosity) 
occur on timescales comparable to the duration of the individ- 
ual pulses (the period of rotation of the neutron star). The pre- 
viously found correlations (based on flux variations on much 
longer timescales) between the cyclotron line energy and the 
X-ray flux are reproduced: that is a negative correlation for the 
super-Eddington transient sources 4U 01 15+63 and V 0332+53, 
and a positive correlation for the sub-Eddington source Her X- 1 . 
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